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We explore models of hard-rod fluids with a finite number of allowed orientations, and construct 
their bulk phase diagrams within Onsager's second virial theory. For a one-component fluid, we show 
that the discretization of the orientations leads to the existence of an artiflcial (almost) perfectly 
aligned nematic phase, which coexists with the (physical) nematic phase if the number of orientations 
is sufBciently large, or with the isotropic phase if the number of orientations is small. Its appearance 
correlates with the accuracy of sampling the nematic orientation distribution within its typical 
opening angle. For a binary mixture this artificial phase also exists, and a much larger number of 
orientations is required to shift it to such high densities that it does not interfere with the physical 
part of the phase diagram. 
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I. INTRODUCTION 

Understanding the phase behavior of colloidal suspen- 
sions of rodlike particles requires an accurate description 
of their microscopic properties. Fluids of hard rods may 
be considered as the simplest systems on which the mod- 
els incorporating particle orientational degrees of freedom 
can be tested fll, l3] . One of the exact theoretical results 
dates back to Onsager ;3] who analyzed the transition 
from a uniform isotropic phase to an oricntationally or- 
dered nematic phase in a fluid of monodisperse hard nee- 
dles. He realized that the average pairwise (rod-rod) ex- 
cluded volume is reduced in the nematic phase compared 
to that in the isotropic phase, and argued that the result- 
ing gain of free volume (and hence translational entropy) 
compensates the loss of orientation entropy (due to the 
nematic ordering) at sufficiently high concentrations of 
rods [3- Onsager derived a nonlinear integral equation 
for the orientation distribution function, a key quantity 
of the theory, which is constant in the isotropic phase 
and peaked about the director in the nematic phase. He 
circumvented the problem of explicitly calculating the ne- 
matic orientation distribution function (ODF) by adopt- 
ing a variational Ansatz, which was numerically checked 
to be rather accurate later . 

The generalization of the Onsager model to binary 
mixtures of rods showed the possibility of strong fraction- 
ation 0, 1^ and even nematic- nematic demixing at suf- 
ficiently high density, driven by a competition between 
orientation entropy and ideal mixing entropy 0, 0, 0- 
The functional forms of ODF's in these studies were ei- 
ther variational Gaussian |^ Q , truncated expansions in 
Legendre polynomials Q, or numerically determined on 
an angular grid 6j or on a scaled angular grid 8]. In all 
these cases the focus was on describing the system with 
a continuum of orientations. 

An alternative approach is to study models with a fi- 
nite number A'^ of allowed orientations while the positions 
of the centers of mass of the rods remain continuous. The 



first such model was proposed by Zwanzig with ori- 
entations of a rod to be restricted to A = 3 mutually 
perpendicular directions = {!,... , A^}- Despite its 
inability to resolve the orientational structure of the one- 
particle distribution function in any detail, it has been 
successfully applied to explore wetting phenomena near 
a single hard wall and in a slit or phase diagrams 
of polydisperse systems 0, The main advantage 

of such discrete models in comparison with the contin- 
uous ones is their computational simplicity. The com- 
bination of spatial inhomogeneity and/or polydispersity 
with a continuum of orientations is rather demanding nu- 
merically, and the computational efforts can be reduced 
significantly by discretizing the orientations pdl . IT^ . 
The hope has, of course, always been that with an in- 
crease of the number of allowed orientations one would 
smoothly approach the continuum limit. Here we show 
that this is not the case. 

The possibility of a continuous interpolation between 
results for the discrete models on the one hand and 
Onsager-like solutions on the other has first been ques- 
tioned by Straley ^| in studies of models with dodecahe- 
dral (A^ — 6) and icosahedral (A^ — 10) symmetries. He 
concluded that they do not trend towards the continuum 
solution due to the single allowed orientation within the 
typical opening angle (~ 7r/9) of the nematic distribu- 
tion at coexistence. Unfortunately, one cannot proceed 
the sequence of models with A^ = 3, 6, 10 any further, 
since a larger fully symmetric set of orientations on the 
unit sphere does not exist. In order to be able to study 
the effect of discretizing the allowed orientations we give 
up part of the symmetry of the set, and this allows us to 
connect continuous and discrete models. We apply our 
method not only to a one-component system of rods, but 
also to binary mixtures, which may be considered as the 
simplest polydisperse systems. 

This paper is organized as follows. In Sec. |n]we de- 
rive the grand potential functional for a model with dis- 
crete number of allowed rod orientations from the On- 
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sager functional. We calculate bulk equations of state 
for specific orientational sets, and determine the num- 
ber of orientations required to resemble the continuous 
Onsager solution. In Sec. IIIII we apply the method to 
construct bulk phase diagrams of binary mixtures of thin 
and thick rods. We demonstrate that their structure can 
be significantly modified due to orientational discretiza- 
tion. A summary and some discussion of our results will 
be presented in Sec. IIVI 



II. MONODISPERSE RODS 

Consider a fluid of hard rods of length L and diameter 
D {L ^ D) in a macroscopic volume V at temperature T 
and chemical potential fx. The "Onsager" grand potential 
functional of the one-particle distribution function 
/9(w) can be written, within second virial approximation, 

as [g 

+ ^ J du}du}'E{u;,Lj')p{u})p{u}'), (1) 

where /3 = {kT)~^ is the inverse temperature, v is the 
rod's thermal volume, and E{uj,u}') is the excluded vol- 
ume of rods with orientations a> and u:' . The function 
p{<2)) is normalized as n = J dCjp{Cj), with n the bulk 
number density (which depends on /3/i). The minimum 
condition 5Vt[p{uj)] / 5 p{u3) = on the functional leads to 
the nonlinear integral equation 



ln[p(w)j/] 
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to be solved for the equilibrium distribution p(w). 

Models with a discrete number N of allowed rod ori- 
entations can be systematically derived from the contin- 
uous model (Q) by dividing the unit sphere into solid 
sectors Aw^, (z — 1, . . . ,N) around vectors (2>j, and fix- 
ing the rod density = p{'^i) within each sector 
as well as the excluded volume i?((2),a>') = E[u3i,Cjj) 
for every pair of sectors. The grand potential func- 
tional ri[pi] of such " orientationally discretized" fluid 
with the density pi ~ p{u}i)Au}i and the excluded vol- 
umes Eij — E[Q}i^Cjj) is 
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with normalization n — X^iLiPi- The last term in Eq. 
© represents the contribution due to the discretization 
procedure into the grand potential fi, i.e. the intrin- 
sic difference between continuous and discrete models. 



For a homogeneous distribution of vectors Cji on the unit 
sphere and AtD^ = Aa> (i.e. for the models with = 3, 6 
and 10 01), it trivially shifts the chemical potential 
Ppd — /?/^ -l- In AtD, which does not have any consequence 
for the solutions pi at a fixed n, and for the thermo- 
dynamics of the isotropic-nematic transition. However, 
when Aa>j is not the same for all i, it acts as an external 
orientational field that tends to favor the larger sectors 
over the smaller ones. This becomes explicit if we con- 
sider the Euler-Lagrange equations that corresponds to 
the discrete functional (or equivalcntly the analog of Eqs. 
©) 
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now to be solved for pi. Note that the equation of state 
p = p{n, T) does not pick up an additional term from 
discretization, 
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but the distributions pi to be inserted into it do depend 
on the discretization. 

Further discussion requires a specification of the set of 
allowed orientations u}i,i — 1, . . . ,N and the associated 
solid sectors At2>i . Unfortunately, it is impossible to com- 
pletely cover a surface of the unit sphere by equal regular 
spherical M-polygons, where M indicates the number of 
polygon's sides (only 5 Platonic solids exist). But sym- 
metries of the function p{u}) can be explicitly included 
into the set of vectors uji in order to simplify the prob- 
lem. For the present study we fix the z axis of the coordi- 
nate system to be parallel to the nematic director n and 
assume uniaxial symmetry of the function p(t2») = p{0), 
with 6 — arccos(a> • n) the angle between a> and n. The 
azimuthal angle is denoted by (f), and hence we charac- 
terize a vector a) — (sin6' cos sin0 sin cos 0) by the 
angles 9 and (j). The "up-down" symmetry of the ne- 
matic phase reduces the orientational space to half the 
upper-hemisphere, i.e. G [0, 7r/2] and £ [0, tt]. As 
we do not expect any azimuthal symmetry breaking we 
restrict attention to A^^ uniformly distributed values for 
(j) for every allowed 9. We have considered a uniform 
distribution of Ng polar angles 9 S [0,7r/2], i.e. 
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as well as a uniform distribution of Nq values of cos 9 G 
[0,1], i.e. 
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Ne-l 
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k^l,...,Ne, l^l,...,N^ (7) 
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with the conventional definition of Aa>i = J^^ sm OdOd4>. 

Figure^shows the dimensionless pressure p* = f3pL^D 
as a function of the dimensionless bulk density n* = 
nL'^D for the grid ® with different Ng and = 5. The 
plateaux (of the sohd fines) correspond to the isotropic- 
nematic coexistence, obtained by equating pressure and 
chemical potential in the two phases. For Ne < 9 the 
transition occurs between the isotropic phase (/) and an 
(almost) perfectly aligned nematic phase (A) with the 
pressure being close to the the ideal gas pressure. Note 
that such grids correspond to a single 9 within the " On- 
sager" opening angle, i.e. < Oi < n/lS < 92- As soon 
as < 01 < 6*2 < 7r/18, or Ng > 9, the distribution 
function p{9) at isotropic-nematic coexistence starts to 
converge to the continuous solution. These results are 
in full agreement with the previous explanation of Stra- 
ley 01 . Equations of state for the model O are very 
similar to Fig. ^but start to resemble Onsager-like dis- 
tribution function for Ng > 80 due to the poor sampling 
of u)i near the nematic director. Equations of state for 
the Zwanzig {N = 3), dodecahedral {N = 6) and icosahe- 
dral {N — 10 ) m odels were calculated using the original 
formulations p^ . and are included for comparison. Our 
results for Ng < 9 seem to converge well to these existing 
results. 

For Ng > 20 the pressure of the high-density nematic 
phase clearly demonstrates a linear dependence on the 
bulk density, i.e. /3p(n) ^ n. With increasing Ng it 
gradually approaches a limiting scaling behavior Pp{n) = 
3n, established for the continuous Onsager solution by 
means of a scaling argument [l6l |. 

The discretization of the rod's allowed orientations 
shows the existence of an "artificial transition" from 
a less-ordered nematic phase (N) to a near-perfectly 
aligned phase (A) for Ng > 9, as indicated by the dashed 
horizontal lines in Fig^ It occurs due to the same com- 
petition between excluded volume and orientational en- 
tropy as in the IN transition, and puts an additional con- 
straint on the description of the nematic bulk state by 
restricted-orientation models. Below we argue that it has 
important consequences for discrete models of (polydis- 
perse) hard-rod mixtures, where separation into nematic 
phases with different composition occurs at sufficiently 
high densities. 



III. BINARY MIXTURES 



+ 2 J dCjdCj'E„^^.{u,-Cj')p,{Cj)p„,{Cj')ii) 

cr.cr' — 1 

with normalization na = f dLbp^jCj). It is known from 
previous work 0, 0, ^| that the bulk phase dia- 
gram of this system exhibits (i) strong fractionation 
at isotropic-nematic (/ — N2) coexistence, (ii) nematic- 
nematic (A^i — N^) coexistence ending in a consolute 
N1N2 point at sufficiently high pressure, and (iii) an 
IN1N2 triple point. The discrete version of this model 
follows directly from ^ as 

2 N 

P^iiPai}] ^YYp<^^{ lii[p^iiy„] - 1 - (3^J.^^ 

(T = l 1 = 1 

^ 2 N 2 N 

+ 2 X! Y E„i.^,jp^ipa'j -^^Pai^^^i^l^u:, (9) 

(T(T' — 1 — 1 (7 — 1 i — 1 

with the densities pai = pcr('^i)^'^ii the excluded vol- 
umes Eci-a-'j = Ecra'{<^i,<^j) and the number densities 
normalization = X^i^i P<yi- Figure [3 shows the phase 
diagrams for discrete systems in the p* — x representa- 
tion with the dimensionless pressure p* = /3pL^Di and 
the mole fraction of thick rods x = ?^2/(^^l + ^2), for 
several orientational grids © with Ng — 11 (a), 20 (b), 
30 (c), 50 (d) and N^ = 10. Note that all four dis- 
cretizations are such that they reproduce the physical 
"Onsager" -like I — N transition at cc = and x = 1, 
at pressures p* ~ 17.7 and 4.4, respectively. However, 
the existence of the artificial aligned nematic phase A 
gives rise to spurious I — A, Ni — A and N2 — A phase 
equilibria, where /, Ni and N2 are the physical isotropic 
and nematic phases. (For the coarsest discretization with 
Ng — 11 (Fig|2Ia)) the A^i phase is stable in a very nar- 
row region beyond the resolution of the picture.) Upon 
refining the discretization from Ng = 11 the IN2A triple 
point (V) shifts to higher pressures, and combines with 
the INiA triple point (A) at Ng = 30 (Fi^c)) to form 
the physical IN1N2 triple point (now denoted by V) and 
an artificial N1N2A triple point (A) at slightly higher 
pressure. Further grid refinements to Ng — 50 yields the 
physical phase diagram with an IN1N2 triple point and 
consolute point as in Ref.|l3], but with a spurious 
N — A coexistence at pressures beyond the "physical" 
part of the phase diagram. At these pressures one does 
not distinguish Ni and N2 nematic phases, but the single 
nematic phase denoted by N here. 



Consider a binary mixture of thin (Di) and thick (-D2) 
hard rods of equal length L and the diameter ratio d = 
D2/D1 = 4 in a macroscopic volume V at temperature 
T and chemical potentials pi and /Z2, respectively. The 
"Onsager" grand potential functional for this system can 
be written as [J| 

2 „ 

/3f][K(u>)}] = dCbp^{Cb)(\n[p,{Lb)i,] - l-(3p„) 



IV. SUMMARY 

We have explored the connections between continuous 
and restricted orientation models of monodisperse and 
binary hard-rod fiuids (in the Onsager "needle" limit 
L ^ D). Our main finding is that a discretization of 
the orientations leads to existence of a nonphysical al- 
most perfectly aligned nematic phase (A) at high densi- 
ties. If the discretization is coarse, i.e., the number of 
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allowed orientations is small, then the A phase can co- 
exist with the isotropic phase (/), and at sufficiently fine 
discretization with the nematic phase N. We also found 
that the continuum limit requires a finer orientation grid 
for a mixture than the one-component fluid. 

Clearly, our results are strongly influenced by the 
adopted limit L ^ D, but we would expect similar ef- 
fects, although weaker, for finite L/D ratio. 

These findings could be relevant for the study of inho- 
mogeneous and/or polydisperse fluids of rods, which are 
computationally more demanding and hence impose the 
use of rather coarse grid of orientations in order to be 
tractable and practical. The present study shows that 
care must be taken with such rather coarse grids, since 



they can give rise to an artificial, discretization-induced 
aligned nematic phase. 
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FIGURE CAPTIONS 

1. Equation of state for models with different num- 
ber of allowed polar angles Nq. Positions of the 
phase transitions are indicated by the horizontal 
lines. The continuous Onsager solution (0ns) can 
be reproduced in the present density interval with 
Ng > 50. Dotted lines correspond to equations of 
state for the Zwanzig (Z), dodecahedral (6S) and 
icosahedral (lOS) models [T^. The dashed horizon- 
tal lines correspond to a spurious nematic- nematic 
(TV — A) transition due to poor discretization of the 
allowed orientations. 

2. Bulk phase diagrams of a binary thin-thick mix- 
ture of hard rods (diameter ratio D2/D1 — 4.0, 
equal length L ^ D2) , in the pressure-composition 
representation, with p* the dimensionless pressure 
and X the mole fraction of the thicker rods. We 
distinguish the low-pressure isotropic phase (/), 
high-pressure nematic phases (A^i and A^2), aligned 
phase A, upper (A) and lower (V) triple phase co- 
existence and an A^i — N2 critical point (*). The 
grey regions, enclosed by the binodals, denote the 
two-phase regimes, and the tie lines that connect 
coexisting phases, are horizontal. 
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